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Neutron stars are extremely relativistic objects which abound in our universe and yet are poorly 
understood, due to the high uncertainty on how matter behaves in the extreme conditions which 
prevail in the stellar core. It has recently been pointed out that the moment of inertia /, the 
Love number A and the spin-induced quadrupole moment Q of an isolated neutron star, are related 
through functions which are practically independent of the equation of state. These surprising uni- 
versal 7 — A — Q relations pave the way for a better understanding of neutron stars, most notably 
via gravitational-wave emission. Gravitational-wave observations will probe highly-dynamical bina- 
ries and it is important to understand whether the universality of the I — X — Q relations survive 
strong-field and finite-size efi'ects. We apply a Post-Newtonian-Affine approach to model tidal de- 
formations in compact binaries and show that the / — A relation depends on the inspiral frequency, 
but is insensitive to the equation of state. We provide a fit for the universal relation, which is 
valid up to a gravitational wave frequency of ~ 900 Hz and accurate to within a few percent. Our 
results strengthen the universality of 7 — A — Q relations, and are relevant for gravitational-wave 
observations with advanced ground-based interferometers. We also discuss the possibility of using 
the Love-compactness relation to measure the neutron-star radius with an uncertainty < 10% from 
gravitational-wave observations. 

PACS numbers: 97.60. Jd 04.30.Db,04.25.Nx, 
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I. Introduction. Neutron stars (NSs) are extremely 
compact objects formed as the end-state of the collapse of 
massive stars, which populate the universe either in iso- 
lation or in binaries. NS-NS binaries are one of the most 
promising sources for second-generation, ground-based 
detectors of gravitational waves, such as Adv. Virgo [T], 
Adv. LIGO j2j and KAGRA |3 (see also the proposed 
third-generation detector ET [3]); in addition, NSs are 
copious radio and X-ray emitters, and they can poten- 
tially be used as laboratories for high-energy and fun- 
damental physics, to probe the behaviour of matter in 
extreme conditions ^5 . A persistent obstacle against ex- 
ploring the full potential of NSs physics lies precisely in 
our ignorance on their inner structure, and in particu- 
lar in the uncertainties of the equation of state (EoS) of 
matter at ultranuclear densities. The star radius R, mass 
M, moment of inertia / and deformability, as measured 
by the tidal Love number A and spin-induced quadrupole 
moment Q, all depend sensitively on the EoS. The lack of 
knowledge on the EoS therefore affects our understand- 
ing of the NS properties and prevents model-independent 
tests of gravity with these objects. 

Despite the multitude of modern EoS proposed in the 
literature, leading to different NS configurations, Yagi 
and Yunes (hereafter, YY) recently discovered some uni- 
versal relations between the moment of inertia, the tidal 
Love number and the spin-induced quadrupole moment. 



which are essentially insensitive to the NS EoS [7]. 
These tantalizing I — X — Q relations open the interesting 
prossibility of breaking degeneracies between these pa- 
rameters and effectively make NSs viable laboratories for 
fundamental physics and astrophysics. 

The tidal Love number used to derive the I — X ~ Q 
relations was computed assuming a static, spherically 
symmetric star placed in a time-independent external 
quadrupolar tidal field. Effectively, this means that these 
relations were derived for stars in isolation. However, 
gravitational-wave detectors in the next years will ob- 
serve coalescing compact binaries, i.e. they will detect 
gravitational waves emitted by NSs at small orbital sep- 
arations, d/R < 10, and high frequencies /gw > 40 Hz; 
thus gravitational-wave observations will probe highly 
dynamical NSs which strongly interact with their com- 
pact companion (either a black hole or another NS). The 
question then arises as to whether the I — X — Q rela- 
tions are actually useful or even correct in situations of 
physical interest. 

The purpose of this Letter is to derive EoS- independent 
/ — A relations which hold true throughout almost the 
entire inspiralling phase of the binary coalescence, and 
to provide accurate fits describing this relation at each 
orbital frequency. 

II. Evaluation of the tidal Love number. Tidal de- 
formability properties of NSs can be described in terms 



of a set of parameters, the Love numbers [HHH], which 
relate the mass multipole moments of the star to the (ex- 
ternal) tidal field multipole moments. In particular, the 
dominant contribution to the stellar deformation is en- 
coded in the electric, I = 2 Love number, which we simply 
call tidal Love number A, and is defined by the relation: 



Qij — AGi- 



(1) 



where Qij is the traceless quadrupole moment of the star, 
and Cij — e?Q^e|^^e?g^ef Ni?a^^5 is the tidal tensor which 
induces the deformation (e? > is a parallely transported 
tetrad, attached to the deformed star). 

Two approaches are currently used to evaluate the 
tidal Love number: a stationary and a dynamical ap- 
proach. In the stationary approach used by YY [5Vll4j. 
the compact bodies forming the binary system are as- 
sumed to be very far apart. Using spacetime perturba- 
tion theory [T5] to study the I = 2 stationary perturba- 
tions of a NS induced by a test tidal field, the quadrupole 
and tidal tensors are evaluated; the Love number is then 
computed from Eq. ^ . As discussed in [SI [H] , this ap- 
proach assumes that the timescale of the stellar defor- 
mation is much smaller than timescales associated to the 
orbital motion, an assumption which becomes less accu- 
rate in the last stages of coalescence. 

In the dynamical approach [TBI Il2|) the evolution of 
the tidal deformation of NSs in compact binaries is mod- 
eled combining the post-Newtonian (PN) description of 
the two-body metric and of the orbital evolution, with 
an affine description of the NS as a deformable ellipsoid, 
subject to its self-gravity, to internal pressure forces and 
to the PN tidal field of the companion. The deformed 
NS is described in terms of five dynamical variables: the 
principal axes of the ellipsoid, and two angles describing 
the orientation of the principal frame; these quantities are 
determined by solving a set of ordinary differential equa- 
tions in time, coupled with the PN equations of motion. 
This approach, called Post-Newtonian Affine (PNA), al- 
lows to compute Qij (t) and Cij (t) in terms of the dynam- 
ical variables, so that the tidal Love number can be eval- 
uated during the inspiral. To parametrize the dynamical 
evolution of the system, it is convenient to use the orbital 
frequency /, instead of time or radial distance. The ratio 
between quadrupole and tidal tensors is then a function 
(the tidal Love function) A(/) [16|,ll7j. and the Love num- 
ber obtained in the stationary approach corresponds to 
the zero-frequency (i.e. infinite orbital separation) limit 
of this function. 

The PNA approach also allows to compute the mo- 
ments of inertia /^ = / x (a^/i?)^, where i = 1,2,3 in- 
dicate the star principal axes (i = 1 corresponds to the 
axis pointing toward the companion), / and R are the 
moments of inertia and the radius of the spherical star, 
and Oi is the z-th axis of the deformed, ellipsoidal star. 
During the inspiral, /i increases, while I2 and I3 decrease. 
///. Results. Using the PNA approach, we have com- 
puted the normalized Love function A = X/M^ and the 



normalized moment of inertia corresponding to the axis 
pointing toward the companion, / = Li/M^, as functions 
of the orbital frequency /; M is the total mass of the 
binary. 

We have performed simulations of NS-NS binaries for 
three different EoS which are expected to cover a wide 
range of NS deformability, APR4, MSI and H4, and masses 
in the range [1.2 -=- 2]Mo . The EoS APR4 describes soft 
NS matter and yields models with high compactness and 
small deformability, whereas MSI describes stiff matter 
and large deformability; H4 provides intermediate config- 
urations. All EoS are modeled by parametrized piecewise 
polytropes as proposed by Read et al. [Ol. 

Our results are summarized in Fig. [I] On the three left 
panels, we plot / versus A for three different values of the 
gravitational wave frequency, /gw = 2/ = 170, 500, 875 
Hz, for equal-mass NS-NS binaries with different EoS. 
The data have been fitted with the following function 

ln/ = ao-|-ailnA + a2(lnA)^ + a3(lnA)^ + a4(lnA)^, (2) 

where the fitting parameters a^ are functions of /gwj 
and are listed in Table [H The dashed lines in the left 
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10-1 9.54 


10-2 -9.72 


10-3 4.12 


10-* 


any 1.95 
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TABLE I. Fitting parameters of the 7 — A relation given by 
Eq. ([2|, for several values of the gravitational wave frequency. 
These fits reproduce our data to within 2%, cf. Fig. [l] The 
last row corresponds to the fit ([3| that reproduces data at 
any frequency to within 5% [cf. Fig. [2] . 

panels of Fig. [T] are the fits corresponding to the selected 
frequencies. 

On the upper, right panel in Fig. [l] the relative error 
{I— I fit)/ 1 fit is plotted versus A, for the selected frequen- 
cies. This error is always < 2%. In the lower panel the 
ratio J(/)//o is plotted versus A, where /q is the asymp- 
totic value of / when the stars are in isolation. This figure 
shows that, as the stars approach the merger, their mo- 
ments of inertia change with respect to the asymptotic 
value, and grow as much as 10%-30%, depending on the 
EoS: for stiffer EoS, the variation with respect to the 
values at infinity is larger. 

Nonetheless, the relative errors [I — If a) /I fa are small 
and only mildly dependent on the EoS, suggesting that a 
simple frequency-independent relation can be found be- 
tween / and A. We find that 
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1.95 - 0.373 In A + 0.155(lnA)2 

- 0.0175(lnA)3 -^0.000775(lnA)^ 



(3) 
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FIG. 1. (Left) The I — X relation is plotted for equal mass NS-NS binaries, three values of the gravitational wave fequency 
/gw = 2/ = (170,500,875) Hz and for the EoS APR4 (x), MSI (O), H4 (D). Markers refer to numerical data, while dashed 
lines are the fits ([2| . (Right Top) Relative fractional errors between fits and numerical results. (Right Bottom) J — A relation 
with moment of inertia normalized by its value at infinity. 



describes very well our numerical results. In the upper 
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FIG. 2. (Top) Fitting curve ([3| (dashed line) and numerical 
results of the 7-A relation, for data-set including points up 
to /gw = 875 Hz and the EOS APR4 (x), MSI (Q), H4 (D). 
(Bottom) Relative fractional errors between fits and numer- 
ical results. Black squares and red circles refer to the fit of 
Eq. ([3|, and to the analytical relation found by YY, respec- 
tively. 

panel of Fig. [2J we compare the fit ^ with our numerical 
results, for the full set of binaries and frequencies we have 
considered. In the lower panel we plot the relative errors 
between the numerical results and the universal fit (l3|, 
showing also how well YY's fit performs in the dynamical 
case. Our fit (Eq. (|3|) reproduces the J— A relation to 
within 5% at any frequency < 875 Hz and for all EoS and 



masses we have considered, while the YY fit becomes less 
accurate as the frequency increases, with fractional errors 
which become of the order of 10%. 

IV. Discussion. NS-NS binaries are the prototypi- 
cal sources for upcoming second-generation gravitational- 
wave detectors. Strong-field and finite-size effects are im- 
portant to model the waveform during the latest stages 
of the inspiral. Our results show that the I — X rela- 
tions discovered by YY in the low frequency regime, can 
be extended to describe the dynamical evolution of NSs 
during the late stages of the inspiral. These results open 
the possibility for strong-field, model-independent tests 
of NS properties. If, for instance, advanced gravitational 
wave detectors LIGO/Virgo measure the tidal Love num- 
ber in a compact binary coalescence to within (5 — 10)%, 
as estimated in [51 [15] , this would allow for an indirect 
estimate of the moment of inertia with roughly the same 
precision. This measurement would be independent of 
(and competitive to) the estimates coming from pulsar- 
timing observations [19] . In addition, as shown in YY, 
these estimates would allow to set constraints on modi- 
fied theories of gravity. 

Although in this paper we have presented only NS-NS 
binaries, our approach also describes the dynamical evo- 
lution of mixed black hole-NS systems as well. We have 
computed the / — A relation for mixed binaries, finding 
similar universal relations during the entire inspiral. 

In this work we have studied the I — X relation, to 
understand the effects of the tidal interaction when the 
stars are at short orbital distance. Spin effects have been 
neglected. They have been considered by YY in the slow 
rotation, low frequency limit. It would be interesting to 
establish whether a simple EoS- independent, universal 
relation exists, between the tidal Love number and the 



spin-induced quadrupole moment Q in the fast rotation, 
high frequency regime. This matter will be investigated 
in a following work. 

We conclude this discussion with some considerations 
on the relation between the tidal Love number and the 
NS compactness C = M/R. YY showed that this rela- 
tion is more EoS-dependent than the I — X — Q relations. 
However, they included in their study hot and young NSs, 
which are unlikely to be members of a coalescing binary 
system. If we consider only old and cold NSs, we find 
that the C — A relation acquires a remarkable universal- 
ity N By computing A in the low frequency limit, for the 
EoS APR4, MSI, H4, and masses in the range [1.2-=- 2] M© 
we find that C is well described by the fit 

C = 3.71 X 10"^ - 3.91 X 10^2 In A -I- 1.056 x 10~^(ln A)^ 

This fit gives the compactness with a relative error < 2% 
The C — A relation can be extremely useful to ex- 
tract information on the NS EoS from a detected grav- 
itational wave signal emitted in a binary coalescence 
If the tidal Love number is extracted by Advanced 
LIGO/Vhgo with an error ainx = (J\/\ ^ 60% [7] 
we can determine the compactness with an error ac '^ 

'^l^^^^^^cJFi^J^^al^^ < 10% C (where we have as- 
sumed a fit < max(C — Cfn)). A much more optimistic 
estimate of the error on the measure of the tidal Love 
number, ainx '^ 5%[TH], would imply a relative error on 
the compacteness of the order of ^ 2% (this remarkable 



decrease of the relative error, can be traced back to the 
~ R^ dependence of A). Since the same detection would 
allow for an accurate estimate of the NS mass, we would 
then know the NS radius with an accuracy of ^ 90% or 
smaller. It should be noted that current estimates of NS 
radius based on astrophysical observations (see [21] and 
references therein), with a claimed error of ~ 10%, are 
highly debated in the literature since they may depend 
on the way the NS surface emission is modeled [35]. A 
measurement of the NS radius, based on gravitational 
wave observations and on the C — X relation, would have 
the same, or better, accuracy and it would be model inde- 
pendent. Such a measurement would be extremely useful 
to put constraints on the NS EoS [23] . 
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^ This is consistent with the results of Ref. 1201 . who studied the 
Q — C relation using a set of EoS describing old, cold NSs, finding 
hints of universality. If we combine the Q — C and the Q — \ 



relation discovered by YY, the universal behaviour of the C — A 
relation naturally emerges. 



